Solving nonlinear multicommodity flow problems by the analytic center cutting plane method

The paper deals with nonlinear multicommodity flow problems with convex costs. A decomposition method is proposed to solve them. The approach applies a potential reduction algorithm to solve the master problem approximately and a column generation technique to define a sequence of primal linear programming problems. Each subproblem consists of finding a minimum cost flow between an origin and a destination node in an uncapacited network. It is thus formulated as a shortest path problem and solved with Dijkstra's d-heap algorithm. An implementation is described that takes full advantage of the supersparsity of the network in the linear algebra operations. Computational results show the efficiency of this approach on well-known nondifferentiable problems and also large scale randomly generated problems (up to 1000 arcs and 5000 commodities

Timing of Millisecond Pulsars in NGC 6752: Evidence for a High Mass-to-Light Ratio in the Cluster Core

Using pulse timing observations we have obtained precise parameters, including positions with about 20 mas accuracy, of five millisecond pulsars in NGC 6752. Three of them, located relatively close to the cluster center, have line-of-sight accelerations larger than the maximum value predicted by the central mass density derived from optical observation, providing dynamical evidence for a central mass-to-light ratio >~ 10, much higher than for any other globular cluster. It is likely that the other two millisecond pulsars have been ejected out of the core to their present locations at 1.4 and 3.3 half-mass radii, respectively, suggesting unusual non-thermal dynamics in the cluster core.Comment: Accepted by ApJ Letter. 5 pages, 2 figures, 1 tabl

Long-term observations of the pulsars in 47 Tucanae - II. Proper motions, accelerations and jerks

This paper is the second in a series where we report the results of the long-term timing of the millisecond pulsars (MSPs) in 47 Tucanae with the Parkes 64-m radio telescope. We obtain improved timing parameters that provide additional information for studies of the cluster dynamics: a) the pulsar proper motions yield an estimate of the proper motion of the cluster as a whole ($\mu_{\alpha}\, = \, 5.00\, \pm \, 0.14\, \rm mas \, yr^{-1}$, $\mu_{\delta}\, = \, -2.84\, \pm \, 0.12\, \rm mas \, yr^{-1}$) and the motion of the pulsars relative to each other. b) We measure the second spin-period derivatives caused by the change of the pulsar line-of-sight accelerations; 47 Tuc H, U and possibly J are being affected by nearby objects. c) For ten binary systems we now measure changes in the orbital period caused by their acceleration in the gravitational field of the cluster. From all these measurements, we derive a cluster distance no smaller than $\sim\,$4.69 kpc and show that the characteristics of these MSPs are very similar to their counterparts in the Galactic disk. We find no evidence in favour of an intermediate mass black hole at the centre of the cluster. Finally, we describe the orbital behaviour of the four "black widow" systems. Two of them, 47 Tuc J and O, exhibit orbital variability similar to that observed in other such systems, while for 47 Tuc I and R the orbits seem to be remarkably stable. It appears, therefore, that not all "black widows" have unpredictable orbital behaviour.Comment: 21 pages in journal format, 9 figures, 4 tables, accepted for publication in MNRAS, several clarifications made and typos fixe

A Shapiro delay detection in the binary system hosting the millisecond pulsar PSR J1910-5959A

PSR J1910-5959A is a binary pulsar with a helium white dwarf companion located about 6 arcmin from the center of the globular cluster NGC6752. Based on 12 years of observations at the Parkes radio telescope, the relativistic Shapiro delay has been detected in this system. We obtain a companion mass Mc = 0.180+/-0.018Msun (1sigma) implying that the pulsar mass lies in the range 1.1Msun <= Mp <= 1.5Msun. We compare our results with previous optical determinations of the companion mass, and examine prospects for using this new measurement for calibrating the mass-radius relation for helium white dwarfs and for investigating their evolution in a pulsar binary system. Finally we examine the set of binary systems hosting a millisecond pulsar and a low mass helium white dwarf for which the mass of both stars has been measured. We confirm that the correlation between the companion mass and the orbital period predicted by Tauris & Savonije reproduces the observed values but find that the predicted Mp - Pb correlation over-estimates the neutron star mass by about 0.5Msun in the orbital period range covered by the observations. Moreover, a few systems do not obey the observed Mp - Pb correlation. We discuss these results in the framework of the mechanisms that inhibit the accretion of matter by a neutron star during its evolution in a low-mass X-ray binary.Comment: 4 figures, 2 tables, accepted for publication in the Astrophysical Journa

Long-term observations of the pulsars in 47 Tucanae. I. A study of four elusive binary systems

For the past couple of decades, the Parkes radio telescope has been regularly observing the millisecond pulsars in 47 Tucanae (47 Tuc). This long-term timing program was designed to address a wide range of scientific issues related to these pulsars and the globular cluster where they are located. In this paper, the first of a series, we address one of these objectives: the characterization of four previously known binary pulsars for which no precise orbital parameters were known, namely 47 Tuc P, V, W and X (pulsars 47 Tuc R and Y are discussed elsewhere). We determined the previously unknown orbital parameters of 47 Tuc V and X and greatly improved those of 47 Tuc P and W. For pulsars W and X we obtained, for the first time, full coherent timing solutions across the whole data span, which allowed a much more detailed characterization of these systems. 47 Tuc W, a well-known tight eclipsing binary pulsar, exhibits a large orbital period variability, as expected for a system of its class. 47 Tuc X turns out to be in a wide, extremely circular, 10.9-day long binary orbit and its position is ~3.8 arcmin away from the cluster center, more than three times the distance of any other pulsar in 47 Tuc. These characteristics make 47 Tuc X a very different object with respect to the other pulsars of the cluster.Comment: Accepted for publication by MNRAS, 18 pages, 11 figure

Population history from the Neolithic to present on the Mediterranean island of Sardinia: an ancient DNA perspective

Recent ancient DNA studies of western Eurasia have revealed a dynamic history of admixture, with evidence for major migrations during the Neolithic and Bronze Age. The population of the Mediterranean island of Sardinia has been notable in these studies –} Neolithic individuals from mainland Europe cluster more closely with Sardinian individuals than with all other present-day Europeans. The current model to explain this result is that Sardinia received an initial influx of Neolithic ancestry and then remained relatively isolated from expansions in the later Neolithic and Bronze Age that took place in continental Europe. To test this model, we generated genome-wide capture data (approximately 1.2 million variants) for 43 ancient Sardinian individuals spanning the Neolithic through the Bronze Age, including individuals from Sardinia{’}s Nuragic culture, which is known for the construction of numerous large stone towers throughout the island. We analyze these new samples in the context of previously generated genome-wide ancient DNA data from 972 ancient individuals across western Eurasia and whole-genome sequence data from approximately 1,500 modern individuals from Sardinia. The ancient Sardinian individuals show a strong affinity to western Mediterranean Neolithic populations and we infer a high degree of genetic continuity on the island from the Neolithic (around fifth millennium BCE) through the Nuragic period (second millennium BCE). In particular, during the Bronze Age in Sardinia, we do not find significant levels of the {“}Steppe{” ancestry that was spreading in many other parts of Europe at that time. We also characterize subsequent genetic influx between the Nuragic period and the present. We detect novel, modest signals of admixture between 1,000 BCE and present-day, from ancestry sources in the eastern and northern Mediterranean. Within Sardinia, we confirm that populations from the more geographically isolated mountainous provinces have experienced elevated levels of genetic drift and that northern and southwestern regions of the island received more gene flow from outside Sardinia. Overall, our genetic analysis sheds new light on the origin of Neolithic settlement on Sardinia, reinforces models of genetic continuity on the island, and provides enhanced power to detect post-Bronze-Age gene flow. Together, these findings offer a refined demographic model for future medical genetic studies in Sardinia

Affine su(2) fusion rules from gerbe 2-isomorphisms

We give a geometric description of the fusion rules of the affine Lie algebra su(2)_k at a positive integer level k in terms of the k-th power of the basic gerbe over the Lie group SU(2). The gerbe can be trivialised over conjugacy classes corresponding to dominant weights of su(2)_k via a 1-isomorphism. The fusion-rule coefficients are related to the existence of a 2-isomorphism between pullbacks of these 1-isomorphisms to a submanifold of SU(2) x SU(2) determined by the corresponding three conjugacy classes. This construction is motivated by its application in the description of junctions of maximally symmetric defect lines in the Wess-Zumino-Witten model.Comment: 41 pages, 1 figure (the published version